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Magneto-optical transitions between Landau levels can provide precise spectroscopic information 
on the electronic structure and excitation spectra of graphene, enabling probes of substrate and 
many-body effects. We calculate the magneto-optical conductivity of large-size graphene flakes using 
a tight-binding approach. Our method allows us to directly compare the magneto-optical response 
of an isolated graphene flake with one aligned on hexagonal boron nitride giving rise to a periodic 
superlattice potential. The substrate interaction induces band gaps away from the Dirac point. In 
the presence of a perpendicular magnetic field Landau-level like structures emerge from these zero- 
held band gaps. The energy dependence of these satellite structures is, however, not easily accessible 
by conventional probes of the density of states by varying the back-gate voltage. Here we propose the 
magneto-optical probing of the superlattice perturbed spectrum. Our simulation includes magneto- 
excitonic effects in hrst-order perturbation theory. Our approach yields a quantitative explanation 
of recently observed Landau-level dependent renormalizations of the Fermi velocity. 

PACS numbers: 73.22.Pr, 71.70.Di, 81.05.ue, 71.70.-d 


I. INTRODUCTION 

The progress in fabrication of well-characterized 
graphene structures led to the availability of gr^hene 
devices with mobilities as high as 6 • 10® cm^/V^. As 
transport properties are strongly influenced by disorder, 
the use of suitable substrates such as hexagonal boron ni¬ 
tride (hBN) are key to substantially reduce the amount 
of substrate-induced bulk disorder compared with con¬ 
ventional Si 02 , dramatically improving the electronic 
and transport properties of graphen^. These substrates, 
however, c an su bstantially modify the physical properties 
of graphen^^. For example, perfectly aligned graphene 
on hBN features a moire pattern of 14 nm periodicity due 
to a small lattice mismatch between the two material^. 
The superlattice potential leads to the opening of a gap 
at the Dirac cone and to the folding of the Brillouin 
zone. Intersections between the original cone and its 
back-folded replica give rise to additional band ga^ en¬ 
ergetically above and below the main Dirac cone^. In 
the presence of a magnetic field perpendicular to the 
graphene, secondary satellite Landau levels emanate near 
these band gaps. Moreover, the large real-space superlat¬ 
tice unit cell allows observing the Hofstadter butterfljEl 
at laboratory accessible field strengtlP*^. 

Recent experimentd^lZiai probing the satellite structures 
and the Hofstadter butterfly of graphene on hBN in a 
magnetic field have employed a varying back gate volt¬ 
age (I4g) to control the effective Fermi level Ep in the 
graphene sheet. Such a measurement does not, how¬ 
ever, provide unambiguous information on the density 
of states: through quantum capacitance effects and de¬ 
viations from the idealized linear density of states of 
graphene (caused, e.g., by the satellite structures), the re¬ 
lation between energy and applied back gate voltage may 
deviate substantially from expectation^. An attractive 
experimental alternative is magneto-spectroscopjff^HUl 


This technique provides information on the energy dif¬ 
ference between the ground state and the particle-hole 
excitation and, thus, between the different Landau lev¬ 
els. Effective Fermi velocities can be extracted from a fit 
of Dirac Landau levels to the appropriate transitions. In 
the presence of a substrate induced moire potential, the 
magneto-optical signal can also provide information on 
the satellites and the evolution of associated Landau lev¬ 
els. The observed shift of the inter-Landau level transi¬ 
tions of graphene on different substrates can be described 
by a renor malized Fermi velocity and attributed to many- 
body effect^iiliSI. Velocity renormalization was a lso m ea- 
sured in magneto-Raman scattering experiment 

In this work, we determine the magneto-optical con¬ 
ductivity of large (140 x 120 nm^) graphene flakes where 
the influence of residual edge effects is suppressed by im¬ 
posing a Berr y-Mo ndragon potential^ at the boundary 
(see also RefsWe evaluate the optical conductiv¬ 
ity by calculating the dipole transition between the dif¬ 
ferent eigenstates of the flake, which we obtain within 
the third-nearest neighbor tight-binding approximation. 
We compare the behavior of a pristine graphene flake 
with a graphene flake aligned with hBN represented by 
a superlattice potentiaP. Absorption lines associated 
with the optical transitions between Landau levels of 
the Dirac fermions show the expected square root de¬ 
pendence on the magnetic held. For graphene on hBN, 
we observe a Hofstadter butterhy on top of each transi¬ 
tion line. The satellite structures above and below the 
Dirac cone, which are the signatures of the Brillouin zone 
folding in graphene aligned with hBN, are shown to con¬ 
tribute to the optical conductivity: our tight-binding ap¬ 
proach predicts additional structures in the optical con¬ 
ductivity that evolve linearly with magnetic held, open¬ 
ing a pathway for optical characterization of graphene- 
substrate interactions. Furthermore, we calculate the 
shift of the inter-Landau levels transition lines of pris- 
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tine graphene in the magnetic field due to an attrac¬ 
tion of excited electron and hole, i.e. magneto-excitons, 
on a tight-binding level. We solve a two-body problem 
and evaluate the direct Coulomb electron-hole interaction 
within first-order perturbation theory. The resulting ve- 
loci ty ren ormalization agrees well with the experimental 
with comp lementary theoretical approaches 
for bulk graphen d^^bsi One advantage of the present 
tight-binding based approach is the inclusion of spatially 
varying long-range potentials due to a substrate. 


II. MAGNETO-OPTICAL RESPONSE OF 
PRISTINE GRAPHENE 

The tight-binding (TB) method is very useful for 
simulating large-scale structures when more sophisti¬ 
cated techniques relying on periodic boundary condi¬ 
tions and small supercell sizes [e.g., density functional 
theory (DFT)] prove computationally challenging due 
to the presence of a magnetic field or additional large- 
scale potential variations. While for small magnetic fields 
(i.e. i? « 0) the density of states (DOS) [see Figj^a)] 
is dominated by size quantization and properties of the 
flake boundaries, for increasing field, Landau levels begin 
to emerge resembling the behavior of bulk graphene. The 
transition from the linear DOS of bulk graphene at B = 0 
to the Landau level regime is governed W the magnetic 
length Ib — 25.5/-\/i3[T][nm]. For flakeLandau levels 
appear when Ib becomes smaller than the flake diameter 
D. The Landau level N for massless Dirac fermions of 
graphene satisfieJ^^ 

= sgn{ N) V F V2ehNB, ( 1 ) 


Eq.§ by optical transitions between Landau levels. Dis¬ 
tortions of the linear DOS by traps or localized charge 
states especially at the device edges, which influence the 
Vbg coupling, are largely eliminated. 

The selection rule for optically allowed inter-Landau 
levels transitions is AN = |A^/| — |A^i| = ±1, where 
Nf{Ni) is the final (initial) Landau level quantum 
numbei!^. Deviations of the realistic graphene band- 
structure from the ideal Dirac cone such as trigonal warp¬ 
ing or the perturbation by the interaction with the sub¬ 
strate will break this selection rule, and transitions with 
|AA^| > 1 will become allowecP^. However, the proba¬ 
bility of these transitions remains small. The main tran¬ 
sition lines visible in the measured infrared absorption 
spectr£iiill 2 l are Ti = — Eq , T 2 = E 2 — E^^ and 

T 3 = Eg — E^ 2 - The analytical forms of the transition 
energies, which correspond to the peak positions in the 
associated absorption spectrum, follow from Eq.Q as 

Ti = ^2\e\hvlB, 

T 2 = yj2\e\nv^pB{V2 + l), (3) 

Tg = ^j2\e\Hvj,B{V3 + V2). 


Employing the eigenstates and eigenenergies of the tight- 
binding calculation for the graphene flake, it is straight¬ 
forward to obtain the optical dipole transitions between 
the different eigenstates due to the coupling to the elec¬ 
tromagnetic field of the infrared laser. The associated 
optical conductivity is given in the length gauge bji^ 




/(eg) - fjeb) 

Ca - eb - + iV ° 


eb)\{a\x\b)\\ (4) 


which coincides with the high density regions in the 
calculated DOS [see black dashed curves in FigliKa)]. 
Here vp is the Fermi velocity of the Dirac dispersion 
Ep = TufIpI. For the present set of tight-binding 
parameterPSl Fermi velocity of pristine graphene is 
v°F = 0.78 • 10® m/s. This value of Fermi velocity deter¬ 
mined from DFT 2 I 1 ^ however, does not include the many- 
body effects seen in experimentd^MSU corrections dis¬ 
cussed below. 

Landau levels in graphene were observed in transport 
experiments and were measured by changing the back 
gate voltage Vbg- The back gate voltage changes the num¬ 
ber of charge carries n = aVbg in graphene and, therefore, 
the Fermi energy 

Ep = hvF^/im ( 2 ) 

through capacitive coupling between the back gate and 
the graphene sheet, where a, the so-called lever arm, is 
a measure for the coupling strength. Combining Eqs. Q 
and ([^ yields straight lines for Landau levels in the Vbg- 
B plane in contrast to the square-root behavior in the 
E-B plane for Landau levels except of iV = 0. By con¬ 
trast, magneto-spectroscopy allows for direct probing of 


where w is a photon frequency of the infrared laser, S is 
the area of the flake, and Cg (e;,) are the eigenenergies of 
the flake corresponding to the eigenstate |o) {\b)). The 
summation indices a and b extend over all the eigenstates 
of the flake in the selected energy window. The matrix 
element (a|i:| 6 ) gives the transition dipole moment be¬ 
tween the two eigenstates. We use the Fermi distribution 
/(eg.h) at zero temperature with the chemical potential 
located at the Dirac point. The difference of Fermi dis¬ 
tributions, /(cg) — /(eh), is non-zero for transitions be¬ 
tween occupied and unoccupied states corresponding to 
particle-hole excitation. Zeros in the denominator signify 
resonant absorption of photons with hw = ea — et- In or¬ 
der to avoid numerical instabilities we have introduced 
a line broadening of 77 = 0.1 meV of the discrete eigen¬ 
states of the flake, which is, however, small compared to 
the physical width of the coarse-grained DOS and, thus, 
of no consequence for the numerical results. 

Real and imaginary parts of the magneto-optical 
conductivity calculated using Eq.Q for the pristine 
graphene flake [FigQb, c)] confirm the dominance of the 
optically allowed transitions of the Dirac cone also for 
the hnite-size graphene flake. The maxima of the (ab¬ 
sorptive) real part of the optical conductivity as well as 
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FIG. 1. (a) The density of states of a 140x120 nm^ pristine graphene flake as a function of energy and magnetic field. 

Black dashed curves denote the Landau levels of massless Dirac fermions (Eq. 0 ). (b) Real (absorptive) and (c) imaginary 
(dispersive) part of the optical conductivity cr (color scale) of the flake (evaluated using Eq. 0) as a function of the excitation 
energy AE = fiw and magnetic field. The black dashed curves are the analytical prediction for the optical inter-Landau levels 
transitions [see Eq. §]■ 


the nodal lines of the (dispersive) imaginary part agree 
well with the analytical prediction (Eq. <§) using the 
Fermi velocity Vp [see black dashed curves in Fig0b, 
c)]. The observed magneto-optical transition lines Ti, T 2 , 
T 3 etc. are narrow with widths of the Landau levels in 
the DOS, of about « 2meV [Fig0a)]. Dipole-forbidden 
transitions with AN > 1 are not visible in Fig0b, c) 
because of their much weaker oscillator strength. 


III. MAGNETO-OPTICAL RESPONSE OF 
GRAPHENE ON hBN 


The tight-binding approximation allows us to study 
the magneto-optical conductivity of more complex struc¬ 
tures with broken periodicity or large periodic super¬ 
cells. Specifically, perfectly aligned graphene on hexago¬ 
nal boron nitride (hBN) features a 14 nm peri odic moire 
pattern due to the small lattice mismatcfP^^. An effec¬ 
tive potential reproducing the main features of the mea¬ 
sured DOS is of the forrrP 


VBN(fO = W{Nj-(Jz + V{r)-\\E lTMB(B,r) • tr^, (5) 


where Uz is the Pauli matrix. This potential has the 
periodicity of the moire pattern. The first term W{f^ ■ az 
is responsible for the breaking of the sublattice symmetry 
and for opening of a band gap of A = 14 meV near the 
Dirac point. The presence of a gap modifies the Dirac 
Landau levels to 


eS{b) 


sgn{N)^2\e\hvj,\N\B + (A/2)2 : N ^ 0 
±A/2 : A = 0±, 

( 6 ) 


The two zeroth Landau levels to which we assign the 
quantum numbers 0 ^ and 0 “ corresponding to energies 
±A/2 emerge as a result of the valley splitting. The val¬ 
ley degeneracy of the other Landau levels is not lifted. 
The second term in Eq. 0 is a smoothly varying poten¬ 
tial accounting for stronger binding in regions where one 
sublattice of graphene is on top of boron and another 
sublattice is in the middle of the BN hexagon. This 
term accounts for the increased electron-hole asymme¬ 
try beyond the third nearest-neighbor tight-binding ap¬ 
proximation. In addition, V (r) may modify the Fermi 
velocity. For the parameters of V(f) considered, the 
Fermi velocity of graphene subject to superlattice po¬ 
tential of hBN decreases to = 0.76 • 10® m/s, which 
is slightly smaller than that of pristine graphene Vp. Fi¬ 
nally, the last term ■ az in Eq. 0 describes, 

on a phenomenological level, magnetic-field dependent 
many-body effects and accounts for the magnetic-field 
induced valley splitting. The lifting of the four-fold de¬ 
generate (two-fold for valley and two-fold for spi n) ze- 
roth Lan dau level was observed in graphene on 
hBlSPSSI and in suspended graphen^^. In graphene on 
hBIvP, the valley splitting was experimentally found to 
be linearly proportional to the magnetic field, a result 
which is still unexplained by theory. Two alternative 
theoretical methods have been previously proposed to 
explain this observation, (i) One model, based on the 
continuous Dirac model in the Hartree-Fock appro xima- 
tion, predicts a Vb scaling of the valley splittin^^MU, 
While this scaling differs from the linear scaling, the re¬ 
sulting numerical values for the valley splitting resem¬ 
bles those of the experiment within the range of inves¬ 
tigated magnetic field strength, (ii) The second model, 
based on the effect of lattice distortion and the inter- 
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action with the substrate, predicts linear scaling of val¬ 
ley splitting with B but underestimates the strength of 
the splitting by an order of magnitud^^. In the present 
work, we account for the valley splitting within the single¬ 
electron model by applying the phenomenological poten¬ 
tial which depends linearly on the magnetic 

field with the spatially averaged value W^mb(^) = 8 i?[T] 
meV/T consistent with experimenlP^. The spatial varia¬ 
tion of the potential Wmb follows the Gaussian shape of 
l/(rO, i.e. WMB{B,f) = Wmb{B) ■ V[^/ (The 
Zeeman spin splitting, also proportional to B, is ne¬ 
glected in the following). 

The DOS of the graphene flake on hBN [Figj^a)] 
as well as the optical conductivity [Figj^b, c)] feature 
the square-root scaling for the N ^ 0 Landau levels 
^ -x/B of massive Dirac fermions given by Eq.([^ 
[see black dashed curves in Figj^a)] for magnetic fields 
up to B « 10 T. By contrast, the linear magnetic field de¬ 
pendence of the zeroth Landau level is determined by the 
effective many-body potential WMB(tO- P Wmb = 0, the 
zeroth Landau levels, i.e. 0"^ and 0“, coincide with the 
magnetic field independent vertical black dashed curves 
in Figj^a)]. The position of the absorption lines in the 
calculated optical conductivity agrees with the energy 
difference between Landau levels (Eq. §) for iV ^ 0 us¬ 
ing Fermi velocity For transitions involving the ze¬ 
roth {N = 0) level, e.g. the Ti transition, we have to 
include in Eq. ([^ the linear magnetic field dependence of 
the 0“ Landau level [see red dashed curve in Fig^b)]. 
The sensitivity of the Ti transition to the many-body val¬ 
ley (and spin) splitting in the zeroth Landau level offers 
an alternative to activation gaps measurements.^ 

At higher magnetic fields, the Hofstadter butterflj® be¬ 
comes visible as ’’diamond”-like structures on top of each 
Landau level in the DOS [see Fig|^a) at B > 10 T]. This 
phenomenon arises due to the competition between two 
length scales: the magnetic length Is = 25.5/-\/B[T][nm[ 
and the scale of the hBN superlattice with a period of 
14 nm. At B = 25 T, i.e. when the magnetic flux 
through the superlattice unit cell is equal to the mag¬ 
netic flux quantum, the DOS features the diamonds with 
the largest extent. The optical conductivity, likewise, 
features a Hofstadter butterfly modifying the absorption 
lines [see Figj^b, c) at B > 10 T]. The energetic width of 
the diamond^WHM) of the Ti transition at 25 T is 20 
meV corresponding to the size of the first Landau level 
diamond obtained from the density of states [Figj^a)[. 

Another superlattice induced effect is the presence of 
satellites above and below the Dirac point near E = 
—0.15 eV and E = 0.14 eV at B = 0 [see arrows in 
Figj^a)], which arise due to the avoided crossings be¬ 
tween the Dirac cone and its replica. The correspond¬ 
ing Landau levels of the satellites evolve linearly with 
magnetic field in the DOS. These satellite structures ap¬ 
pear also in the optical conductivity as an absorption 
line starting at AB = 0.3 eV at B = 0 and linearly in¬ 
creasing with B [marked by arrows in Figj^b,c)[. Its 
slope corresponds to the slope of the satellite Landau 


levels. The pronounced difference in the magnetic field 
evolution of Landau levels of the satellites (~ B) and the 
main Dirac cone (~ \/B) is not obvious when measured 
as a function of Vbg. Optical transition spectroscopy 
appears as an attractive alternative for disentangling 
Dirac-like and Schrodinger-like dispersions in the band- 
structure. Magneto-optical experiments for graphene on 
hB]>fH]have, up to now, not focused on the satellite struc¬ 
tures induced by the moire pattern. While challenging, 
the experimental observations of the satellites and their 
magnetic field evolution would yield valuable insight into 
the influence of the moire potential on the bandstructure. 


IV. MAGNETO-EXCITONIC SHIFT OF THE 
OPTICAL INTER-LANDAU LEVELS 
ABSORPTION LINES 


In the single-particle picture we used so far, the cal¬ 
culated optical conductivity yields the position of ab¬ 
sorption lines Ti in agreement with analytical predic¬ 
tions [Eq.([^] based on the difference between the cor¬ 
responding Landau levels of Dirac fermions [Eq.(^ and 
Eq.(|^] In the experimenfl ^^b2 ( however, one observes 
deviations which can be conveniently parameterized in 
terms of transition-line dependent Fermi velocities in 

the AB ^ ^/l3 relation (Eq.([^). Their origin are ob¬ 
viously excitonic interactions correcting for particle-hole 
excitations absent in addition spectroscopy when varying 
the back-gate potential t4g- 

The formation of magneto-exciton d^^ l ^^ l by optical ex¬ 
citations between the Landau levels is a many-body ef¬ 
fect. To lowest order the excitonic wave function is a 
product of an electron and a hole wave function: 

^/?M(^el,ni) = V’Ar(f^el)V’M(fT), (7) 

where fei (Fh) is an electron (hole) coordinate, ipN and 
ipM correspond to the wave functions of the and 
Landau levels undergoing optical transitions. The 
corresponding Dirac equation reads 

Hi) = (Fel, Fi), (8) 


where 


B = ^ VF^ ■ {pi - eA) - 


i—el,h 


47reoe [Fei - rh| ’ 


(9) 


e is a dielectric constant and eg is the permittivity of the 
vacuum. We use the value of e = 5 found for graphene on 
Si 02 within the random phase approximation (RPA[P^ 
and also for graphene on hBlNp^. In the absence of the 
Coulomb interaction [second term in B, Eq.([^] the exci¬ 
tonic energy coincides with Eq.([^, i.e. the single-particle 
energy difference. The Coulomb interaction changes the 
transition energies and leads to the observed energy shift 
due to particle-hole attraction. Magneto-excitonic ef¬ 
fects were studied in detail for the two-dimensional elec- 
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FIG. 2. Same as Figj^ but for a 140x120 nm^ graphene flake deposited on hexagonal boron nitride (hBN) giving rise to a 
periodic snperlattice potential with period 14 nm period. The black dashed cnrves in (a) correspond to the Landan levels of 
Dirac fermions with a finite mass [Eq. §]. The red dashed lines indicate the magnetic field dependent valley splitting in the 
zeroth Landan level due to the additional phenomenological potential WMB{B,f). The satellites above and below the Dirac 
point at B = 0 and their magnetic held evolution are marked by arrows. In (b) and (c) black dashed curves denote the optical 
transitions Ti, T 2 and T 3 , predicted by the difference between the corresponding Landau levels for massive Dirac fermions 
(Eq.®). The red dashed trace in (b) corresponds to the Ti transition when the magnetic held dependent valley splitting of the 
zeroth Landau level is included. The optical transitions between Landau levels at B = 0 of satellites and their B-held evolution 
are marked by arrows (compare to Fig®b,c)). 


tron ga^sHSH (2DEG) and for graphen^i^EsIlSl pg^j._ 
ticular, it was shown that the excitation energy con¬ 
sists of several contributions: (i) the (non-interacting) 
single-particle exciton energy AE = — E^, (Eq.(|^); 

(ii) the direct Coulomb two-particle interaction between 
the particle and the hole; (iii) the annihilation and cre¬ 
ation of electron-hole pairs at different points of the Bril- 
louin zone; and (iv) the exchange interaction. The direct 
Coulomb term (ii) is negative and gives rise to the exci- 
tonic binding. Therefore, the optical transition energy is 
reduced relative to the estimate in the single-particle pic¬ 
ture. Contributions (iii) and (iv) provide positive higher- 
order corrections, slightly reducing the Coulomb attrac¬ 
tion. Since this effect is already, to a certain extent, 
empirically accounted for by the dielectric response of 
the medium e, we neglect an explicit treatment of these 
terms. 


Within the tight-binding approximation applied to 
graphene dots, it is possible to treat magneto-excitonic 
effects. For comparison with the experiment we include 
the dielectric environment for graphene on a Si 02 sub¬ 
strate by taking into account the dielectric constant e = 5 
when evaluating the Coulomb interaction term. We treat 
the direct Coulomb electron-hole interaction in first-order 
perturbation theory. Consider an optical transition be¬ 
tween the two Landau levels N and M. Each of the Lan¬ 
dau levels comprises a large number of eigenstates of the 
flake [see Fig.® a))]. A photoexcited state can be writ¬ 
ten as a superposition of flake eigenstates with energies 


in a vicinity of the optically excited Landau levels: 

^Nliirex, Eh) = E (10) 

a,b 

The expansion coefficients Cab can be obtained by insert¬ 
ing Eq.([Io| into the Dirac equation Eq. (®: 

(e^ - ef)Cab - Y. yfb'Ca'b' = eCab, (11) 

a' ,b' 


with (Cfe^) being eigenenergies of the flake in the vicin¬ 
ity of (M‘^) Landau level. The Coulomb matrix 
element is 


wab _ 
rn>h> — 


iireoe 

dfeidrh 




( 12 ) 


Tel - rh\ 


Eq. ini) can be viewed as the Bethe-Salpeter equation 
(BSEp3 adapted for a finite-size graphene flake. The 
Coulomb integral Eq.(12) gives rise to the exciton bind¬ 
ing energy. The full solution of Eq. is numerically 
challenging as the size of the subspace coupled by the 
Coulomb matrix is given by the product of the number 
of particle and hole states for each Landau level. 

For each Landau level involved in the transition, we 
consider the energetically closest Ns eigenstates of the 
flake. Civen the size of our model system, we find that 
using Ns « 30—40 (depending on the magnetic field and 
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FIG. 3. (a, b, c) Real part of the optical conductivity [renormalized to its maximum] of the transition lines Ti, T 2 and T 3 as 
a function of photon energy. Conductance traces at different magnetic fields are shifted vertically. Black traces correspond to 
single-particle transitions, green traces correspond to optical transitions taking into account the direct Coulomb electron-hole 
interactions from first-order perturbation theory, (d) The position (dots) and width (bars) of spectral lines extracted from the 
first and second moment of the single-particle (black dots) and excitonic excitation spectra (green dots). The black dashed 
curves correspond to the analytically predicted Landau level transitions given by Eq.O with a Fermi velocity Vp = 0.78 • 10® 
m/s. The green dashed curves correspond to the fit of the lines with Eq. , however, with renormalized Fermi velocities; 

= 0.69 • 10® m/s, Vp^ = 0.745 • 10® m/s and Vp^ = 0.755 • 10® m/s. (e) The ratio between Fermi velocities of T 2 (right) and 
T 3 (left) transitions to the Fermi velocity of Ti: theory versus experiment. 


the Landau level index) is a good compromise between 
accuracy and numerical effort: the matrix dimensions of 
are given by Ng x 7V| « 10®, were the evaluation 
of each matrix element requires a double spatial integral 
with a Coulomb kernel. To keep the problem numeri¬ 
cally manageable, we (i) restrict ourselves to first-order 
perturbation theorj®, i.e. we include only the diagonal 
Coulomb matrix elements in Eq.(ll I and (ii) reduce 
the size of the graphene flake to 24 x 24 nm^, to speed up 
the calculation of each matrix element. As we verified nu¬ 
merically, the diagonal elements indeed, dominate. 
For each electron-hole pair \a,b) the corresponding first- 
order energy correction to the noninteracting electron- 
hole transition energies = (ca — e&) is given by the 
diagonal matrix element 


= vf = — 

ta,b ''ab 


dfeidrh 




Tel - ?’h| 


(13) 


The binding energy of the [a, b) exciton is in the diagonal 
approximation given by ea,b = 

The optical conductivity of each electron-hole pair, 
which is related to the dipole matrix element between 
the electron ]a) and hole \b) eigenstates, is calculated 
using Eq. however, now with the excitonic energy 
contribution included. First-order perturbation theory 
results in a pronounced shift of the Ti, T 2 and T 3 tran¬ 
sition lines in the direction of smaller photon energies, 
when comparing the noninteracting single particle pic¬ 
ture [black traces in Figj^a,b,(^ with the two-particle 
correction [green traces in Fig|^a,b,c)]. The position 
and the width of the spectral lines is determined by their 
first and second moments. Neglecting the excitonic ef¬ 
fect the transition energies [black dots in Figj^d)] agree 
well with the analytical prediction for inter-Landau lev¬ 
els OTtical transitions given by Eq. © [black curves in 
Figj^d)] with a Fermi velocity Vp = 0.78 • 10® m/s de¬ 
termined by our choice of tight-binding parameters for 
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pristine graphene. In contrast, the two-particle excitonic 
corrections shift the lines to lower transition energies Ai? 
[green data points in Fig|^d)]. A fit of these transition 
lines to the analytic prediction AA ~ \/l3 yields now line- 
specific Fermi velocities. In particular, for Ti the Fermi 
velocity is Vp^ = 0.69 • 10® m/s; for T 2 , = 0.745 • 10® 

m/s and for T 3 , Vp^ = 0.755 • 10® m/s. Such a ve- 
locity renormalization, also seen in the experimenlfAE^, 
clearly reflects excitonic effects. To distinguish effects of 
the Landau-level specific velocity renormalization from 
other bandstructure effects that uniformly affect all tran¬ 
sitions, we consider the ratio of the renormalized Fermi 
velocities of T 2,3 transitions to that of Ti. W e find rea¬ 
sonable agreement with the experiment j ^ ^ F 2 J measuring 
optical inter-Landau levels transitions for graphene on 
Si 02 and also with alternative theoretical approaches for 
bulk graphen j^^ l ®® ! [Fig^e)]. The present tight binding 
ansatz allows for including long-range disorder, such as 
puddles observed in graphene on Si 02 . We account for 
the effects of charge puddles by applying a smooth dis¬ 
order potential with an amplitude of 25 meV (50 meV) 
and a correlation length of 10 nm (5nm[P^. Disorder only 
slightly modifies the theoretically calculated line shifts 
and velocity renormalization. The resulting optical inter- 
Landau levels transition lines remain well-defined in the 
presence of disorder in agreement with magneto-optical 
experiments for graphene on Si 02 . 

However, the magneto-optical measurements of 
graphene on hBJ^A^ notably differ from other measure¬ 
ments and theoretical predictions [see right most point 
in Figl e)]. This mismatch arises from a substantially 
(« 30%) lower experimental value of v]} than theoreti¬ 
cal predictions and was attributed to many-body effects 
influencing the zeroth Landau levefAl. Corrections to 
the dielectric response of the material (i.e., the effective 
e) should affect all transitions similarly, and thus can¬ 
not explain the observed large ratio Vp /vp . One effect 
large enough to explain these findings is the substan¬ 
tial splitting of the four-fold degenerate Landau levels 
of pristine graphene, as found for graphene on hBN in, 
e.g., quantum capacitance measurement^. However, any 
such splitting should shift the Ti line to higher photon 
energies [as seen in Figj^b,c)], i.e. in the opposite direc¬ 
tion to that observed in the experiment. Clearly, future 
magneto-optical experiments for graphene on hBN are 


called for to shed more light on this puzzle. 

One remark should be added regarding the large values 
of the magnetic fields used in the simulations to map out 
Landau levels. Since we consider a flake of a smaller size 
than in the previous section, the Landau levels emerge at 
higher magnetic fields. In particular, the Landau level N 
is formed as soon as 2'</2NIb is smaller than the smallest 
dimension of the flake. For example, the N = 3 Landau 
level participating in the T 3 optical transition becomes 
distinguishable only at H > 30 T for the small flake used 
to calculate magneto-excitonic corrections. Clearly, the 
optical response of the larger flakes used in experiment 
can still be extracted from the present calculation by ex¬ 
trapolating the square root behavior of the optical tran¬ 
sitions with transition-dependent Fermi velocity to lower 
fields. 


V. CONCLUSION 

We have simulated the optical properties of graphene 
flakes with and without moire potential for aligned 
graphene on hBN using the tight-binding approximation. 
Our simulations show that the magneto-optical response 
allows probing the satellites due to Brillouin zone fold¬ 
ing and, unlike probing by back gate voltage, to clearly 
distinguish energy levels that scale linearly with B from 
those that feature a square-root scaling, ^ ^/B. We have 
also shown that excitonic effects can be included in a 
tight-binding description. We validate our predictions 
for a Landau-level specific renormalization of the Fermi 
velocity by comparing to experimental data for graphene 
on Si 02 , opening a pathway towards the description of 
excitonic effects in larger structures. For graphene on 
hexagonal boron nitride, current experimental data does 
not fully agree with theoretical predictions, calling for 
further experimental and theoretical studies. 
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